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Abstract 



Even within the framework of the leading logarithmic approximation the eigenvalues of 
the BKP kernel for states of more than three reggeized gluons are unknown in general, 
contrary to the planar limit case where the problem becomes integrable. We consider 
a 4-gluon kernel for a finite number of colors and define some simple toy models for 
the configuration space dynamics, which are directly solvable with group theoretical 
methods. Then we study the dependence of the spectrum of these models with respect 
\ to the number of colors and make comparisons with the large limit case. 

in 
in 

)0 1 Introduction 

o 

In quantum field theory and statistical mechanics the 1/N (or large N) expansion pQ is 
^ | a well known and extensively used perturbative framework whenever the theories under 
investigation present an internal symmetry typically related to groups like SO(N) or SU (N). 

Quantum Chromodynamics is one of the theories mostly studied under this approxima- 
tion even if, as a physical gauge theory, it is characterized by a gauge group SU(N C ) where 
the number of colors N c is just 3. Recently, thanks to the renewed interest induced by the 
ADS/CFT correspondence, the N = 4 SYM theory in the infinity color (planar) limit has 
been intensively studied and several important results achieved. 

The fact that the planar N = 4 SYM is expected by the theoretical community to be 
solvable and that it is dual to a superstring sigma model has led several theorists to look for 
hints, in the absence of any supersymmetry, for the existence of a possible dynamical system 
dual to planar QCD sharing with it some integrability properties. The starting points are the 
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integrable structures unveiled many years ago at one loop in standard perturbation theory 
and some hints of possible integrability at two loops in the planar limit. 

The first evidence of an integrable structure at one loop in QCD was found [2j by L.N. 
Lipatov in the framework of the Regge limit of scattering amplitudes whose behavior may be 
conveniently described by systems of interacting reggeized gluons, as we shall briefly review 
in the next section. The integrable dynamics, associated to the evolution in rapidity of such 
a system, appears when one is taking the large N c approximation, which makes the BKP 
kernel [3] to resemble the structure of an Heisenberg XXX spin chain, but for a non compact 
SX(2,C) "spin". 

Going beyond the large N c approximation, even in the lowest orders in perturbation 
theory in the coupling constant, is a formidable task and it is very difficult also to try to 
estimate the error one faces when computing quantities for infinite N c (planar limit) instead 
of at iV c = 3. 

It is the purpose of this work to introduce some finite toy models, which share the same 
color structure of the BKP systems and can be studied to determine the dependence of 
the spectrum on the number of colors N c . They are characterized by a configuration space 
which is no more the transverse plane but a finite vector space associated to irreducible 
representations of the SU(2) group so that one may use group theoretical methods to analyse 
some of these models. 

This is of course not providing any concrete answer for the question related to the real 
QCD problem, but nevertheless can be of some help. Moreover some toy models may be 
interesting by themselves as dynamical systems. 

We start in the next section with a short review of the properties of the system of 
interacting reggeized gluons in the Leading Logarithmic Approximation. In section three we 
consider the color structure for the four reggeized gluon system and describe how to use a 
convenient basis for it. In section four we construct some finite toy models which are studied 
in some details in a couple of subsections. After the conclusions in few appendices we give 
more details on symmetry structures and on the features of these toy models. 

2 BKP Kernels 

Let us start by giving a brief overview of the kernels which encode the evolution in rapidity 
of systems of interacting reggeized gluons in the leading logarithmic approximation (LLA). 
The reggeized gluons provide a convenient perturbative description of part of the QCD 
degrees of freedom in the Regge limit (also known as the small x limit) and appeared in 
the investigations of the leading dependence of the total cross sections on the center of mass 
energy in the LLA, which is associated to the so called BFKL (perturbative) pomeron [I]. In 
the simplest form, the BFKL pomeron turns out to be a composite state of two interacting 
reggeized gluons "living" in the transverse configuration plane in the colorless configuration. 
The construction of the kernel reflects the property that in the Regge limit the scattering 
amplitude factorizes in the impact factors which determine the coupling of the external 
particles to the t-channel reggeized gluons and in a Green's function which exponentiates the 
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kernel and contains the rapidity dependence of their composite state. Such a dependence 
can be analyzed in terms of the spectral properties of the kernel and in particular one is 
interested in the eigenvalues and eigenstates associated to the leading behavior. Because of 
this the spectral problem is often formulated in quantum mechanical terms with the kernel 
being the "Hamiltonian" and its eigenvalues the "energies" . 

In the case of a colorless exchange the Hamiltonian is infrared finite and in LLA is 
constructed summing the perturbative contributions of different Feynman diagrams: in par- 
ticular the virtual ones (reggeized gluon trajectories) ui and the real ones (associated to an 
effective real gluon emission vertex) V . One writes formally H = u>i + u 2 + T{T 2 V\ 2 where 
Ti are the generators of the color group in adjoint representation . In the colorless case one 
has T{F 2 = —N c and finally one obtains: 

F 12 = ln |pi| 2 + In | P2 | 2 + ^i-ln \p 12 \ 2 Pl p* 2 + -L In \p 12 \ 2 p* lP2 - 4tf(l) , (2.1) 

where \l/(x) = d In T(x)/dx, a factor a s = a s N c /ir has been omitted and the gluon holomor- 
phic momenta and coordinates have been introduced. 

The gauge invariance gives the freedom to choose a description within the Mobius space 
|6] , wherein the functions describing the positions of the two reggeized gluons in the transverse 
plane are zero in the coincidence limit. In this space the BFKL hamiltonian has the property 
of the holomorphic separability (-H12 = ^12 + ^-12)- Moreover a remarkable property is its 
invariance under the Mobius group, whose generators for the holomorphic sector in the 
Mobius space for the principal series of unitary representations are given by: 

M r 3 = Pr d r , M+ = d r , M~ = -^<9 r . (2.2) 

The associated Casimir operator for two gluons is 

M 2 = \M\ 2 = -p 2 12 d 1 d 2 , (2.3) 

where M = Y^ r =x ^ r anc ^ M r = M~ , M^). Due to this symmetry the holomorphic and 
antiholomorphic parts of the Hamiltonian can be written explicitly in terms of the Casimir 
operator: indeed one has, after defining formally J(J — 1) = M 2 , 
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V»(J)+^(1- J) -2^(1). (2.4) 



Labelled by the conformal weights h = + iv , h = + iu, where n is the conformal 
spin and d — 1 — 2iv is the anomalous dimension of the operator O/j ^(p ) describing the 
compound state [?], the eigenstates and eigenvalues of the full hamiltonian in eq. (12. lft . 
H i2E h ~ h = 2x h E h - h , are respectively given by: 



(2.5) 
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and 

(\ + 1 71 1 \ fl + 1 71 1 \ , . , 

Xh = x(", n)=4> [—^ + ™ J + i> ( — 2^ - iv ) ~ 2 ^ ■ ( 2 ' 6 ) 

The leading eigenvalue, at the point n — v = 0, has a value Xmax = 4 In 2 « 2.77259, 
responsible for the rise of the total cross section as s° isXmax , which corresponds to a strong 
violation of unitarity. 

Let us now consider the evolution in rapidity of composite states of more than 2 reggeized 
gluons [3]. The BKP Hamiltonian in LLA acting on a colorless state can be written in terms 
of the BFKL pomeron Hamiltonian and has the form (see [2]) 

H n = —— ^2 T k TiH k i . (2.7) 

c l<k<l<n 

This Hamiltonian is conformal invariant but cannot be solved in general. Nevertheless the 
case of three reggeized gluons, where the color structure trivially factorizes, is solvable [2] 
and different families of solutions were found [ESS]. Physically these states are associated 
to the so called odderon exchange [10] and in particular the family of solutions given in [9] 
corresponds to eigenvalues up to zero (intercept up to one) and are the leading one in the high 
energy limit. Moreover they have a non null coupling to photon- meson impact factors [11J. 

The case of more than three reggeized gluon is in general not solvable but if one considers 
the color cylindrical topology when taking the large N c limit the resulting Hamiltonian 

ttn=\ [#12 + #23 + • • • + H nl \ =h n + h n (2.8) 

is integrable, i.e. there exists a set of other n — 1 operators q r , which commute with it and 
are in involution. They are given, in coordinate representation, by 

ii<i2<— Al- 
together with similar relations for the antiholomorphic sector. In particular, q 2 = M 2 is the 
Casimir of the Mobius group. This is the first case where integrability was found within the 
context of gauge theories analyzing the Green's function in some kinematical limit. This 
integrable model is a non compact generalization of the Heisenberg XXX spin chains (2J [12] 
and has been intensively studied with different techniques in the last decade [13], HH [15] fT6] 

nans]. 

Here we remind the result for the highest eigenvalue of a system of four reggeized gluons 
in the planar, one cylinder topology (1CT), case: H^ip^ = 2i? 4 1CT -?/>4. The maximum value 
found, for zero conformal spin, is 

£ 4 1CT = 0.67416 . (2.10) 

In general for an arbitrary number n of reggeized gluon in the cylindrical topology the 
leading eigenvalues have been found to be positive for even n and negative for odd n and 
asymptotically behaving as 1/n [To! I 
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The following are two important question that are unfortunately very hard to answer: 
what are the eigenvalues at finite N c = 3 and what is in general their dependence in iV c ? 
One may be tempted to apply variational or perturbative techniques to the spectral problem, 
which nevertheless appears to be quite involved. In any case a first step consists of analyzing 
the color structure, which simplifies a bit in the case of four reggeized gluons in a total 
colorless state. 



3 Color structure 

We consider the BKP kernel H 4 for four gluons, given in eq. (12. 7p . This is an operator acting 
on 4-gluon states, which may be represented as functions of the transverse plane coordinates 
and of the gluon colors v({p i }) aia ' 2a:ia " 1 . Let us concentrate here on the color space. 

It is convenient, due to the fact that the four gluons are in a total color singlet state, to 
write the color vector v a i a 2 a 3 a 4 [ n terms of the color state of a two gluon subchannel. Let us 
therefore start from the resolution of unity for a state of two SU(N C ) particles in terms of 
the projectors P{Ri]a{al onto irreducible representations: 

1 = P x + P 8A + Pss + Ao+io + Pvr + Po = Yl W > ( 3 - U ) 

i 

where TrP[Rj\ = di is the dimension of the corresponding representation. Let us note 
that we have chosen to consider a unique subspace for the direct sum of the two spaces 
corresponding to 10 and 10 representations. This is convenient for our purposes and we shall 
therefore consider just 6 different projectors to span the color space of two gluons. 

If we consider gluons (1, 2) to be the reference channel we introduce as the base for the 
color vector space the set {-Pt-RiJaiail} °f projectors and write 



^aia2Q3i4 



5> l ( P ^E) OT « = J3« < P la [i2 i ]. (3.12) 



Note that one could have also chosen other reference channels corresponding to a description 
in terms of projection onto irreducible representations of other gluon subsystems. Having 
chosen a color basis, we find that the next step is to write the BKP kernel with respect to 
it. We can slightly simplify the expression for the kernel since for a colorless state we have 
'Yl li T i v = which implies that T{T 2 v = T 3 T 4 t> (an similarly for the other permutations of 
the indices). Therefore one may write: 
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T\T 2 (H 12 + H u ) + T\T3 (H 13 + H 24 ) + T\T 4 (H u + H 23 ) . (3.13) 



Let us now write explicitly the action of the color operators T{Tj = Yla'P^Pj' which are 
associated to the interaction between the gluons labelled i and j. We start from the simple 
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"diagonal channel" for which we have relation T(Tj = — Y2k a kPij[Rk] with coefficients = 
(N c , ^j-, —1,1). Consequently we can write in the (1, 2) reference base 

(f&vY = -ajV j = - (Av) j , (3.14) 

where A = diag(ak)- The action on v of the TiT 3 and XiT 4 operators is less trivial and 
is constructed in terms of the 6j symbols of the adjoint representation of SU(N C ) group. 
We shall give few details in the appendix A and write directly the results, in terms of the 
symmetric (after a similarity transformation) matrix operators: 

(fAv) 1 = fe * * ?) vi = -(PACvf (3.15) 

and 

(fif^y = ~J2 fe SjCiakCfsA v* = - (SCACS v) j . (3.16) 

The matrix C is the crossing matrix built on the 6j symbols and S = diag(sj) is constructed 
on the parities Sj = ±1 of the different representations Rj. 

We can therefore write the general BKP kernel for a four gluon state, given in eq. f)3.13p . 

as 

Hi = ±r [A (H 12 + H u ) + CAC (H 13 + H 2i ) + SCACS {H lA + H 23 )} (3.17) 

One can check that if we make trivial the transverse space dynamics, replacing the Hij 
operators by a unit operators, the general BKP kernel in eq. (12. 7p becomes H n = f 1 and 
indeed one can verify that A + CAC + SCACS = N c l. 

Let us make few considerations on the large N c limit approximation. As we have already 
discussed, in the Regge limit one faces the factorization of an amplitude in impact factors and 
a Green's function which exponentiates the kernel. The topologies resulting from the large -/V c 
limit depend on the impact factor structure. In particular one expects the realization of two 
cases: the one and two cylinder topologies. The former corresponds to the case, well studied, 
of the integrable kernel, Heisenberg XXX spin chain-like. It is encoded in the relation: 
T{Tj — > which leads to H± = \ {Hn + H 23 + H u + H41). It is characterized by 

eigenvalues corresponding to an intercept less then a pomeron. The latter case instead is 
expected to have a leading intercept, corresponding to an energy dependence given by two 
pomeron exchange. Consequently one expects at finite iV c a contribution with an energy 
dependence even stronger. In the two cylinder topology the color structure is associated to 
two singlets {5 aia2 8 aiai , together with the other two possible permutations). Such a structure 
is indeed present in the analysis, within the framework of extended generalized LLA, of 
unitarity corrections to the BFKL pomeron exchange [Tj5] and diffractive dissociation in 
DIS [20J, where the perturbative triple pomeron vertex (see also [21]) was discovered and 
shown to couple exactly to the four gluon BKP kernel. 
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It is therefore of great importance to understand how much the picture derived in the 
planar N c = oo case is far from the real situation with N c = 3. One clearly expects for 
example that the first corrections to the eigenvalues of the BKP kernel are proportional to 
but what is unknown is the multiplicative coefficient as well as the higher order terms. 



4 Toy models 

In this section we shall consider a family of models, different from the BKP system, which 
neverthelss share several features with it and can be used to judge how the large N c ap- 
proximation might be more or less satisfactory. Moreover these systems may be considered 
interesting by themselves as quantum dynamical systems. 

A state of n reggeized gluons undergoing the BKP evolution, described by the kernel in 
eq. (12. 7p . belongs to a vector space of functions on a domain given by the tensor product of 
the color space 8 n and the configuration space M. 2n , associated to the position or momenta 
in the transverse plane, of the n gluons. Indeed the BKP kernel is built as a sum of product 
of color (T k Ti) and of configuration (H kl ) operators; the latter, on the Mobius space, can be 
written in terms of the Casimir of the Mobius group, i.e. in terms of the scalar product of 
the generators of the non compact spin group SX(2, C): B. ki = H k i(M k ■ Mi). 

We are, therefore, led to consider a class of toy models where the BKP configuration 
space M 2n is substituted by the space V™ where V s is the finite space spanned by spin states 
belonging to the irreducible representation of SU(2) with spin s. In particular we shall 
consider quantum systems with an Hamiltonian fitting the following structure: 



Hn = -W S TkTifiSkA), (4.18) 



c KkKn 



where Si are the elements of the su(2) algebra associated to the particle i in any chosen 
representation and / is a generic function. A particular toy model is therefore specified by 
giving the spin s of each particle ("gluons") and the function /. In the following we shall 
consider two specific cases for the 4 particle system: 

a) A spin s — 1 case in a global singlet state v SiV = 0). If / is the identity map 
than the "spin" configuration dynamics is very similar to the one of the color sector. In 
order to have a system which behaves similarly to the BKP case we first put a constraint on 
the two particle operators, which describe the basic interaction. In particular we consider 
the family of functions 



f a (x) = 2Re 



1 



V> ( g + V -"( 4 + 2x ) 



2^(1). (4.19) 



Remembering that for conformal spin n = the BFKL Hamiltonian is given by H k \ = 
ip I | + y j + (Mfc + Mi) 2 ] — 2-^(1), one immediately recognizes that the f a is associ- 



2Re 

ated to the substitution ~ + L\- — > —a Sf • which assures to have the same leading eigenvalue 
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for any a, since both expressions have the value zero as upper bound. The parameter a 
will be chosen in order to constrain the full 4-particle Hamiltonian (14.181) to have the same 
leading eigenvalue as the QCD BKP system in the large N c limit (at zero conformal spin). 
In this system, the BKP toy model, we shall investigate finite N c effects. 

b) A system TOYAdj.Fund with / the identity function and spin s — 1/2. Such a system 
in the large N c limit in the case of one cylinder topology becomes the well known Heisen- 
berg XXX spin chain system which is integrable. We shall perform some check on the N c 
dependence again for the 4-particle case. 

c) Moreover in order to have another check of the approach we shall also consider a 
model where the 4 particle belong to the fundamental representation of SU(2) for both the 
"color" and the "spin" so that we can perform a comparison with standard results from the 
spectroscopy of isospin-spin systems. We place these checks in the appendix C. 



4.1 BKP toy model 

In order to explicitly study this finite system, described by the Hamiltonian in eq. (14.181) act- 
ing on vector states with dimension (8x3) 4 and singlet under both SU(3)c and SU (2) spin con f, 
it is convenient to choose the color decomposition in 2-particle subchannel irreducible repre- 
sentations described in section 3 and adopt a similar approach also on the "spin" degrees of 
freedom. After that one is left with the problem of diagonalizing an Hamiltonian which is a 
matrix 18 x 18, a problem addressable with any computer. Without the singlet restriction 
on the spin part the problem in general is much more complicated to be easily solved and 
may be addressed in future investigations. 

Let us therefore proceed by introducing for 2 particle spin 1 states the resolution of unity 
1 = Qi + Q3 + Q5 = J2i Q[Ri] which let us write SiSj — — J2k bkQij[Rk] with b k = (2, 1, -1) 
(c.f. with at: first, second and second last terms). It is, therefore, straightforward to write 
from a power series representation (Qij[Rk\ are projectors): 

= X^(-Wi;[#fc] • (4-20) 

k 

Using the corresponding crossing matrices D and the parity matrix S' one obtains relations 
very similar to the one reported in eqs. ( I3.14l) -( l3.16l) . which read 

= f(-b J y = (B V y , (4.21) 

(/ (f&) v)' = {j2Dif(-h)Dlj v* = (DBDvY (4.22) 
and 

(/ (f^) v) 3 = fe s ' 3 D lf(-bk) D's'A v l = {S'DB DS> vf . (4.23) 
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From the above results for the two particle representation basis, we can write the explicit 
form of the Hamiltonian for this toy model, going beyond the one given in eq. (I3.17p . Indeed 
we obtain 

H 4a = ^-(A®B + CAC® DB D + SCA CS ® S'DB DS') (4.24) 

which contains a dependence on N c and on the parameter a through the function f a given 
in eq. 1Q51) . 

Let us note that in the large N c limit one faces for the Hamiltonian two possible cases: 
the one cylinder topology (1CT) which corresponds to the simpler Hamiltonian 



n 



ICT 
4a 



B + S'DBDS' 



(4.25) 



and the two cylinder topology (2CT) corresponding to the even simpler Hamiltonian 



q_/2CT 



Nr. 



N c f(S 1 S 2 )-N c f(S 3 S. 



2B 



(4.26) 



Let us remark that while in the case of N c > 3 we consider a basis for the vector states made 
of P[Ri]Q[Rj] with 18 elements since in the color sector there is also the Pq projector, the 
case N c = 3 is characterized by a basis of 15 elements. 

As already anticipated, in order to study a toy model resembling the spectrum of the 
BKP system of 4 gluons, we require that, in the large N c limit in the one cylinder topology, 
the leading eigenvalue must be the same as the one found for the corresponding integrable 
BKP system, whose value was given in eq. (12.101) . This fact fixes the value of the parameter 
a = 2.80665. We are therefore left with an Hamiltonian which is just a function of the 
number of colors N c . 

Let us now consider its spectrum for the cases N c — 3 and N c = oo. Here we report the 
values followed by their multiplicities. Note than for N c = 3 there are 15 eigenvalues while 
they are 18 for any other value of N c . For the case Nq = oo we specify also the topology 
they belong to. 



/ N c = 3 \ 



7.04193 


(xl) 








5.51899 


(x2) 




( 


N c 


= oo \ 


1.12269 


(x2) 






5.54518 


(x3) 2CT 


-3.89328 


(x2) 






0.67416 


(x3) 1CT 


-4.04744 


(xl) 






-4.27838 


(x3) ICT 


-4.27838 


(xl) 






-7.81242 


(x3) 2CT 


-7.81242 


(xl) 






-8.67983 


(x3) ICT 


-9.18576 


(x2) 




\ 


-10.0168 


(x3) 2CT J 


-12.6743 


(x2) 








-14.1005 


(xl) J 
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We track the flow from Nq — 3 to N c = oo: the first three highest eigenvalues (in bold) 
are moving to the same leading value (in bold) which corresponds to two BFKL pomeron 
exchange (in two cylinder topology). The fourth and fifth highest eigenvalues (underlined) 
are instead moving to the leading eigenvalues of the one cylinder topology case (which are 
three instead of two because of the larger basis for N c > 3). With very good approximation 
one finds that the N c dependence of the leading eigenvalue E is given by 

/ 2 465 \ 

E (N C ) = E (oo) \l + —J . (4.27) 

One can see that for this toy model the large N c approximation corresponds to an error of 
about 27%, an error which is not negligible because the coefficient of the leading correction 
to the asymptotic value, proportional to is a large number. 

It is also easy to investigate the color-configuration space mixing which is encoded in the 
eigenvectors. We report some results in the appendix B. 



4.2 TOY Adj 

.Fund 

We now move to study the toy model described at point (b) at the beginning of section 4, 
again to see how the large N c approximation works. It is described by the Hamiltonian 

T-tAdo,Fund = -jj- ^2 ^ fe ^YY' (4.28) 

c l<k<l<n 

acting on spin singlet states. Again we consider the large N c limit. The one cylinder topology 
is associated to the well known Heisenberg XXX spin chain with Hamiltonian 



1 n 

o/lCT _ 1 ST^ a i (a oq\ 

n -Adj,Fund ~ ^ 2 2 ' V^-OV) 
i=l 



which we shall now consider for the case of n = 4 particle. In this case at large N c we have, 
as before, also the two cylinder topology associated to the Hamiltonian 



ft 2 Adf,Fund = 2 y y • ( 4 -30) 



The spectrum for the one cylinder topology case is well known from Bethe Ansatz meth- 
ods [22] and for total zero spin of a 4-particle spin chain the possible eigenvalues are and 
— 1 (see table II in [23] for J = —1/2 in their notation). The two cylinder topology is 
characterized by the eigenvalues +1/2 and —3/2. 

At finite iV c we rewrite the Hamiltonian in a similar way to the BKP toy model case (see 
eq. (14.241) where the B and D matrices are defined for / the identity map and for the group 
SU(2) in fundamental representation). At N c = 3 it corresponds to a 10 x 10 matrix while 
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10 



15 



21) 



25 



Figure 1: N c dependence of the eigenvalues of the model TOYAdj,Fund- 



for N c > 3 it is given by a 12 x 12 matrix. The leading eigenvalue as function of N c can be 
easily computed 



and indeed goes to the value 1/2 in the large N c limit. Let us note that if one considers the 
planar approximation (in the 2CT configuration), the leading eigenvalue would be underes- 
timated with a relative error of (E (3) — E Q (oo)) / E (3) ~ 40% w.r.t. the case N c = 3. In 
Fig. [I] we report the N c dependence of all the eigenvalues in the range 3 < N c < 25. 

Similar models, but in a higher spin representation, can be constructed in order to main- 
tain the integrability in the large N c limit. One simply needs to consider for any irreducible 
representation s of the particles the function / to be a corresponding specific polynomial as 
described in [22] . 

5 Conclusions 

We have introduced a family of dynamical models describing interacting particles with color 
and spin degrees of freedom. The main motivation was to study within this framework how 
much the large N c approximation is significant when one is trying to extract the spectrum 
of these quantum systems. 

Indeed in some relevant physical cases the only results available are restricted to the case 
with a planar structure resulting from the large N c approximation, when integrability arises 
and gives the possibility to exactly solve the problem. These facts are seen when considering 




(4.31) 



4iV c 
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QCD scattering amplitudes in the Regge limit and LLA approximation, characterized by the 
BKP dynamics. 

We have focused our study to the the case of four particles and considered in details three 
toy models. One toy model (case (c) in section 4.2) was considered to test our computational 
method based on group theory since one is able to make a direct comparison with results 
already known from other methods used in spectroscopy. 

The first model presented in section 4.2 (a) is aimed to mimic to some extend the be- 
havior of the 4 gluon BKP kernel, since we have forced it to have in the large N c limit the 
same leading eigenvalues of the BKP system for both one and two cylinder topologies. We 
were able to compute the different eigenvalues of this toy model as function of N c and we 
have found that the leading one at N c = 3 present corrections of almost 30% w.r.t. the 
planar approximation, which one may understand in terms of a large coefficient in the 1/N% 
correction term. The mixing in color-spin configuration structure has been also studied. 

Another model (case (b) in section 4.2) was considered since in the large N c limit it gives 
rise to the one cylinder topology Heisenberg XXX spin chain which is integrable. For the 
spin 1/2 case we have found at finite N c = 3 corrections to the leading eigenvalue of about 
40%. 

Let us note that in our analysis we were restricting ourselves to study the toy model 
Hamiltonians on the space of states which are singlet with respect to the 577(2) "spin" 
configurations. This was a choice dictated by technical reasons but one should look forward 
to extending the investigation to all the possible states. 

These kind of models and possibly more general ones appear to be interesting also by 
themselves and we feel that they deserve more studies in order to see, for example, if some 
remnant from integrability can be traced back at finite N c . 
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Appendix A 

In this appendix we note a few facts about the crossing matrices introduced in section 3 and 
4 for the SU (N c ) group. Related considerations may be found in [2U ESI EHl EH EH] where 
explicit expressions for the crossing matrices can be found and therefore will not be given 
here. 

Let us rewrite in graphic notation the operator T{Fj in the basis and the color vector 
state in the basis (1,2). 
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= - E* a k 



3 4 

v = E> i>' ] [i 



i J i J 12 

Let us compute the first non trivial crossing case, TiT 3 v, remembering to rewrite the 
final result again in the basis (1, 2). In a graphical notation we have 



1 3 





3 4 






1 3 


1 2 








)o:> 







2 

3 4 



1 



1 2 



1 2 

where the crossing matrix (essentially 6j symbols) can be written as 



c* = 




In a similar way one can also trace the action of the T X T 4 operator. One can see that in 
the last relation there is an asymmetry due to the fact that one divides by the dimension 
of the /c-representation. It is convenient to perform a similarity transformation to work 
with a symmetric crossing matrix. For this purpose it is sufficient to introduce the matrix 
A = diag(dj) and define the new symmetric matrix C — > A^^CA^ which acts on the vectors 



with components v 1 



A' 



i 



Appendix B 



Let us consider the BKP toy model described and analysed in section 4.1. From numerical 
investigations one finds that the leading eigenvector t> and the two closest subleading t>i j2 
at iV c = 3 have the following components 



/ 0.590 


PiQi 


\ 


0.085 


PiQb 




0.344 


PsaQ3 




0.199 


PssQi 




0.199 


P8SQ5 




0.293 


P10+T0Q3 




0.179 


P27QI 




\ 0.574 


PnQ'o 


) 



0.166 
0.342 
0.317 
0.385 
0.267 
0.598 
0.421 



PiQz 

PsaQi 
PsaQ 5 
PssQs 
Pio+ibQi 

P27Q3 



\ 



V2 



( 



\ 



-0.775 


PxQi 


\ 


0.002 


PiQ 5 




0.008 


P%aQ3 




0.123 


PssQi 




0.114 


P&sQb 




0.268 


-P10+10Q3 




0.151 


P27Q1 




0.525 


P27Q5 


J 



As one can see the eigenvector vq of the highest eigenvalue is even which the two fold 
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degerate next lager eigenvalue has eigenstates of both parities {v\ odd and v 2 even). 

In the large N c limit case the eigenvectors of the three fold degenerate leading eigenvalue 
of the 2 cylinder topology are 



/ 



w, 



(2CT) 



( 1 P1Q1 ) *>i 



(2CT) 



1 

3 

1 

f 



-Plo+ioQi \ 
P10+T0Q5 
P27Q3 

P0Q3 ) 



w 



(2CT) 



1 
1 

VE_ 

3^2 
373 

\ 3^ 



Pi 



10+10 
P27Q1 
P27Q5 
P0Q1 
P0Q5 



Q 3 \ 



J 



Again also in this system we can track the same parity properties, which are invariant 
under the flow in N c . 

Similarly one may investigate the states associated to the one cylinder topology at N c = 
00 and their corresponding partners at finite N c . For brevity we just report here the two 
most relevant states in the N c infinity limit: 



w, 



(1CT) 



PsaQi 

PssQ3 
P8AQ5 



w 



(1CT) 



PssQi 
PsaQz 
PssQb 



w. 



(1CT) 



0.245 
0.663 



(P0Q1 
(PoQ 5 



P27Q1) 
P27Q5) 



where Z\ ~ 0.815, z 3 ~ 0.405 and z 5 ~ 0.415. We stress that w^ CT ^ has no corrispective at 
N r = 3. 



Appendix C 

This appendix is devoted to check in one specific case that our approach gives result in 
agreement with other methods widely used in spectroscopy. We start by considering a 
system of n particles in the bifundamental representation of SU(N C ) x SU(2), characterized 
by an Hamiltonian (14. 18ft (with / the identity function) 

c \<k<l<n 

which can be written in terms of the quadratic Casimir operators of SU(N C ), SU(2) and 
SU(2N C ) D SU(N C ) x SU{2) (see [29]). 

Indeed the tensor products of TkSi are among the generators of SU(2N C ), so it is useful 
to introduce the entire algebras for this group 



A; = 1,2,3 = / 
V2T k St k = AT 2 + 3,...,4iV c 2 -l; Z = 1,2,3 



a *=lviri T * k = A,...,N 2 c +2; a = 1, . . . , iV c 2 - 1 (5.33) 
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SU(2) <g> SU(2) 


SU(A) 


(/ii,/i 2 ,/i 3 ) = K S U{±) 


1® 1 


1,20,35 


(0,0,0)(0,2,0)(4,0,0) 


1® 3 


15,45 


(1,0,1)(2,1,0) 


1® 5 


20 


(0,2,0) 


3® 3 


15,20,35,45 


(1,0,1)(0,2,0)(4,0,0)(2,1,0) 


3® 5 


45 


(2,1,0) 



Table 1: Correspondence between irreps of 577(4) and 577(2) ® SU{2) 

: 1/2 5kk'- The Hamiltonian for this system can be rewrit- 



with the normalization Tria^ct^ 
ten as 



7~Laii- 



fund 



AN r . 



1 N 2 
C 2 -2n 



2N r 



s(s + l) 



(5.34) 



where the quadratic Casimir operators C n are defined as in [2H] and s = 1/2. Note that all 
the operators introduced above depend on the irreducible representation of the symmetry 
group to which they refer to. 

We are interested in the real representations so we set N c = 2 and consider the case of 
only four particles. The symmetry group of the model becomes SU{2) ® SU{2) C SU(4) 
and eq. (I5.34p . written for the four particle in a global singlet state, takes the form 



'Hau- 



fund 



cm - ~ 



(5.35) 



In order to find its spectrum the next step consists of analyzing the irreducible representation 
content of each symmetry group of the model. So, for four particle with spin 1/2, one has 
(we specify also the multiplicity) 



2®2<g>2®2 = 2(1) + 3(3) + 5, 

and in the SU(4) case 

4®4®4®4 = 1 + 3(15) + 2(20) + 35 + 3(45) 



(5.36) 



(5.37) 



Then we need to study the SU (2) <g> SU(2) content of these SU (4) irrep. This can be done 
using the results of [30J and in particular the entries of table 1, where the values in the third 
column are Dynkin indices. 

So for particles in a total singlet state (1 eg) 1) the Hamiltonian in eq. (I5.35f) admits 
four eigenvalues, each for a different irrep of SU(4), with a 2- fold degenerate eigenvalue 
corresponded to 20su(4)) (see eq. (15.371) ): 



+ 



15 

16' 

_3_ 

16 

_9_ 

16' 



for irrep 35 
(2x), for irrep 20 
for irrep 1 



(5.38) 



15 



and these are in perfect agreement with the spectrum evaluated with the method used previ- 
ously throughout the paper (we take advantage from the formulas of [HI] for the eigenvalues 
of a quadratic Casimir operator as functions of the Dynkin indices). 

As a final remark we want to emphasize that the method of writing the Hamiltonian in 
terms of the Casimirs can be applied to systems with any number of particles (at the price, 
increasing their number, of a growing complexity in the induced irreducible representations) 
and moreover the analysis may not be restricted to singlet subspaces. Unfortunately it is not 
clear how to define a method for interacting particle not in the bifundamental representation. 
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